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1.0  INTRODUCTION 


The  absorption  of  infrared  radiation  by  water  vapor  in  the  atmospnere  is 

of  great  interest  since  it  plays  a  prominent  role  in  determining  atmospneric 

transmission  to  solar  or  laser  radiation,  and  the  heat  balance  of  the  lower 

atmosphere."  Of  particular  concern  for  long  path  transmission  of  laser  signals 

is  the  so-called  continuum  absorption  in  the  window  spectral  regions.  Tnis 

2 

structurel ess  absorotion  has  been  studied  experimentally  by  Bignell  and  by 
Burch^  and  collaborators;  it  nas  been  treated  phenomenologically  by  Clough, 
et  al^.  For  self-broadening  of  H^O  in  the  1000  cm  ^  window  region,  the  ab¬ 
sorption  is  characterized  by  a  strong  negative  temperature  dependence. 

Theoretical  models  advanced  to  explain  tne  absorption  include  far-wing 
absomtion  from  distant  strong  lines,  dimer  (or  cluster)  absorption,  colli¬ 
sion-induced  absorption,  or  some  combination  of  the  above. 

The  far-wing  theory  and  calculations  presented  in  this  Report  contain  the 
following  features: 

•  The  primary  approximation  is  the  si ngl e-perturber  aporoximation .  This 
assumes  that  the  far-wing  absorption  can  be  obtained  by  calculating  the 
absorption  due  to  a  single  radiator-perturber  pair,  and  then  multiplying 
by  the  total  number  of  pairs.  This  is  consistent  with  the  observed  den¬ 
sity  dependence  of  the  absorption  coefficient.  The  single-perturber  ap¬ 
proximation  for  absorption  in  window  regions  seems  to  have  first  been 
discussed  by  Baranger.^  The  approximation  can  also  be  obtained  as  a  1 i m- 

o 

iting  case  from  the  dipole  autocorrelation  function  in  the  time  domain. 

One  of  the  present  Authors  has  also  discussed  this  approach  using  a  T-ma- 
g 

trix  formalism  directly  in  the  freouency  domain. 
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•  The  present  treatment  considers  absorption  arising  only  from  unboune 
states  of  the  two-body  system. 

•  Line  coupling  contributions,  althougn  formally  contained  in  the  theory, 
have  not  been  included  in  the  present  calculations. 

•  It  is  assumed  that  the  anisotropic  interaction  can  be  treated  using  a 
multipole  expansion.  For  strong,  close  collisions  this  assumption  may 
break,  down. 

•  Perturbation  theory  in  the  anisotropic  potential  is  also  assumed.  Ar 
isotropic,  Lennard-Jones ,  potential  is  included  and  is  treated  exactly. 

•  A  unitarity  cutoff  (analogous  to  Anderson's  well-known  cutoff)  is  em¬ 
ployed  in  the  strong  collision  limit.  A  rigorous  justification  of  this 

cutoff  procedure  is  lacking  and  remains  a  problem  in  terms  of  interpreta¬ 

tion  of  results. 

•  The  theory  rigorously  satisfies  the  fluctuation-dissipation  theorem 

8*10 

(FDT)  if  the  unitarity  cutoff  is  ignored.  If  such  a  cutoff  is  ap¬ 
plied,  the  validity  of  the  FDT  depends  on  the  manner  in  which  the  cutoff 

procedure  is  employed.  In  the  past,  failure  to  satisfy  the  FDT  has  lead 

to  ambigious  i nterpretations  of  far-wing  absorption. 


In  Sections  2,  3,  and  4  we  present  the  general  theoretical  approach  used 
in  the  calculations  of  far-wing  absorption.  In  Section  5  we  present  results 
from  these  calculations. 
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2.0  GENERAL  THEORY  -  THE  SINGLE-PERTURBER  APPROXIMATION 


In  the  problem  of  line  broadening  in  olasmas,  Baranger  has  discussed  ao- 
olication  of  the  “one-electron"  approximation  to  the  far-wing  ana  window  spec¬ 
tral  regions.  For  the  case  of  molecular  collisional  broadening  we  shall  em¬ 
ploy  the  term  "single-perturber"  approximation. 

The  justification  given  by  Baranger  for  the  si ngl e-Derturber  approxima¬ 
tion  ’s  essentially  that  the  far-wing  absorption  will  be  governed  by  the  be¬ 
havior  o*  the  dipole  autocorrel ati on  function  at  very  short  times.  This 
statement  must  be  regarded  with  some  caution  because  the  approximation  is  as 
much  an  expansion  in  density  as  in  small  times. 

The  general  statement  of  the  approximation  is  as  follows: 

% 

"At  sufficiently  low  densities  one  may  calculate  the  far-wing  spectrum  by 
computing  the  absorption  due  to  a  single  pair  consisting  of  one  radiator  and 
perturber,  and  tnen  multiplying  by  the  total  number  of  pairs." 

It  snould  be  noted  that  this  is  distinct  from  the  binary  collision  ap¬ 
proximation  which  assumes  the  absorDtion  can  be  computed  by  calculating  the 
absorption  arising  from  a  single  radiator  and  Npgr  perturbers  which  it  may  in¬ 
teract,  and  then  multiplying  by  the  number  of  radiating  molecules,  The 
binary  collision  approximation  envisions  that  the  radiator  interacts  with  one 
perturber  at  a  time,  however,  over  the  time  scale  of  interest  for  the  radia¬ 
tive  process,  the  radiator  interacts  with  many  perturbers.  At  low  densities, 
the  binary-collision  approximation  leads  to  the  well-known  lorentzian  (impact) 
1 ’reshape  for  isolated  lines  near  the  linecenter.  As  pointed  out  by  Baran- 
ger,^  the  rounding  of  the  Lorentzian  at  the  line  center  is  always  a  many-body 


effect.  The  single-perturber  approximation  is  actually  a  special  case  of  tne 
binary  collision  approximation  in  which  multiple  scattering  can  be  ignorec 
over  the  time  scale  of  interest. 

Baranger  has  also  noted  that,  altnough  the  single-perturber  approximation 
is  not  valid  at  the  line  center  (it  contains  singularities),  there  is  gener¬ 
ally  a  range  of  frequencies  for  which  the  impact  and  single-perturber  approxi¬ 
mations  overlap.  Schematically,  for  an  isolated  line,  the  single  perturbed 
lineshape  function  assumes  tne  form 

i  r  (o)  -  ojfj) 

(I) 

where  is  the  line  center  frequency.  Equation  (1)  is  also  the  approximate 
form  of  a  Lorentzian  several  halfwidths  away  from  resonance,  i.e.  for 
(u-tof1)»r,  except,  in  the  case  of  a  simple  Lorentzian,  the  halfwidth  func¬ 
tion,  FCw-w^),  is  usually  replaced  by  its  value  at  line  center,  r(0).  To  the 
extent  that  the  halfwidth  function  varies  slowly  over  a  frequency  range  of 
several  halfwidths,  there  will  clearly  be  some  overlap  of  the  two  approxima¬ 
tions. 

Our  mathematical  discussion  of  the  single-perturber  approximation  will 
proceed  via  the  frequency  domain  analysis  of  Ref.  9.  In  Appendix  A  we  indi¬ 
cate  how  these  same  results  can  be  obtained  from  a  small-time,  low-density 
limit  of  the  dipole  autocorrelation  function  in  the  time-domain. 


-om  Section  2  of  Ref.  9  the  absorption  coefficient  (cm  )  can  be 


outec  as 


4t2o) 

-  -gg-  X”(W) 


wnore  x"(w)  ca^  be  wn-ppen  ( cons' stent  w*tn  the  PDT) 


X"(w)  »  tanh 


14>M  +  0(-oj)] 


*  (1  -  e~Jho})<t>(u)  .  (3b) 

“  f<^(o>)  -  <3S(-co)J  ,  (3c) 

w'th  S=  1  / k0T .  anc  woene  $(w)  ’s  the  Fourier  transform  of  the  dipole  autocorre- 

8  1 0 

1  at  ion  function  $(t).  (The  definitions  of  x11  in  Refs.  B  ani  9  differ  by  a 
factor  of  n  tne  radiaton  density.) 

An  exact  many-body  formula  *or  $(w)  is 

»  4-  E  e($,)!< liMlF>l*6(^F-^,-hco)  (A) 

u  IF 

Here  2  denotes  the  system  volume,  I,  F  refer  to  exact  eigenstates  of  the 


many-body  system,  and 


<?<«,>  - 


is  the  equilibrium  canonical  density  matrix. 


Now  in  the  si ngl e-perturber  approximation,  we  replace  Eq .  (4)  by 


<t>(u)  =  pr  E  e(E,)|<HM!F>i25(EF-ErM 
u  IF 
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where  I>,  i F>  denote  exact  eigenstates  of  a  two-bocy  Hamiltonian  consisting 
of  a  single  radiator  and  perturDer.  Also  in  Eq.  (6).  N  1  s  the  number  o' 
pairs  given  for  foreign  broadening  by 

Npairs  *  Nrad  Nper  ■  ( 7 ) 


(8) 

(9) 


The  eigenstates  ,  I > ,  F>  are  not  simply  product  states  consisting  or  a 
radiator  state,  a  perturber  state,  and  a  state  describing  the  relative  trans¬ 
lational  motion.  That  is,  in  general  the  internal  states  of  the  two  molecules 
are  coupled,  true  bound  states  of  the  pair  are  admissible,  and  u  is  the  dipole 
moment  for  the  composite  two-body  system. 


In  tne  analysis  which  follows,  we  carry  out  a  perturbative  treatment  in 
the  anisotropic  interaction,  considering  absorptive  contributions  associated 
only  with  unoound  states  of  the  system.  In  this  case,  the  basis  states,  can 
be  taken  as  product  states  (including  an  isotropic  interaction  which  is 
treated  exactly),  and  for  y,  we  can  write,  for  foreign  broadening, 


i*  m /*, rad,  (10) 

and  for  self-broadening 

X  -  Frad  *  '‘per  .  <n> 
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The  difference  between  Eq .  (10)  and  Eq .  (11)  essentially  makes  up  for  the  fac¬ 
tor  of  2  difference  between  Eq .  (7)  and  Eq.  (9)  and  has  the  consequence  that 
se1 f-oroadening  can  be  treated  in  oasically  tne  same  fashion  as  foreign  broa- 
deni  ng  . 


In  carrying  out  the  analysis,  we  write  the  two-body  Hamiltonian  as 


+  V. 


(12) 


where  any  Durely  isotropic  interaction  is  contained  in  Hg,  and  where  V  denotes 
tne  anisotroDic  interaction  between  the  two  molecules.  Here,  anisotropic 
means  that  V  is  ang  ■  e-depenaent  and  couples  the  internal  states  of  the  two 
molecules  with  the  relative  translational  vector  connecting  the  molecules. 


Tne  eigenstates  of  Hg  are  product  states  of  the  form 

H0I'>  *  €j|i>  . 

|i>  ■  1  jjmj >  | JjMj >  | kf- > 


(13) 

(14) 


wne*-e  j.  denotes  a  radiator  state.  J.  denotes  a  perturber  state,  and  k .  de¬ 
notes  a  state  of  relative  translational  motion;  the  m's  and  M's  are  magnetic 
quantum  numbers  fo1"  the  interna1  states.  The  state  k..>  is  an  eigenstate  of 
the  Hami 1  ton i an 


where  r 


’  s  the 


restive  coordinate,  r= 


|ki>  -  ^ 

?  ,  V„(r) 


2  2 

and  k  /2m,  with  m  the  reduced  two-body  mass, 

the  energy  e  ..  can  be  written  as  the  sum  e.=e-  ■*■£ 
density  matrix  factors  according  to 


kj>  (15) 

is  the  isotropic  potential, 
For  the  basis  states  above, 
.  +e,  ,  and  the  unperturbed 

j .  k . 
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©<€j)  =  a(Cji)  e(€jj)  G<€ki>  (16) 

The  T-matrix  formulation  in  Ref.  9  rigorously  incluaes  line  coupling  con¬ 
tributions.  However,  our  initial  calculations  nave  ignorec  this  complication 
In  this  regard,  we  also  remark  that  some  of  the  equations  presented  in  Section 
2.4  of  Ref.  9  are  actually  a  slight  oversimpl i f ication  of  the  uncoupled  line 
approximation.  The  essence  of  the  uncoupled  line  approximation  is  tnat  when 
one  encounters  a  matrix  element  product  of  the  form 

<  jj^j  |  Jl\  jf  mf  ><  if  ,mf ,  |  £1  jj  ,m } ,  >  (17) 

one  sets  j'=j.  and  however  one  does  not  set  m..  ,=m^  or  m*,=mf  for  tne 

magnetic  quantum  numbers. 

For  $(w) .  we  obtain 

* 

<t>{u)  -  r  E  +  1)j  <jj|  | ft J  |jf>  j*  *  *jf(w  -  ojfj)  .  (18) 

In  the  above  equation: 

1.  The  prefactor  r  is  the  same  renormal  ization  factor  discussed  in  Eq. 
(2.13)  of  Ref.  8,  and  in  Ref.  9,  i.e. 

r  -  Tr(e  ^H°)/ Tr(e”^H)  .  (19) 

In  all  of  our  calculations,  we  have  approximated  this  factor  by  unity. 

2.  The  reduced  matrix  element  is  defined  by  the  Wigner-Eckhart^  theorem 
in  the  form 

<jjmii^|jfmf>  -  <ji||fi||jf> 

•  E  em(iflmfmlif1ijmi)  (20) 

m«0,±l 

The  reduced  matrix  element  satisfies  the  symmetry  relation 
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(21) 
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From  £q.  (22)  we  can  note  that 

iZi) 

wnere  fiw..=(er-E  ■  )•  By  letting  «--<*>  in  Eq .  (IS),  interchanging  names  of  dummy 
summations  and  Dy  making  use  of  symmetry  oroperties  (21).  (27).  and 

(28),  one  may  readily  show 

<>(-od)  =  <t>(u)  (29) 

Equations  (27)  ana  (29)  are  the  funaamental  relations  whicn  guarantee  that  tne 
FDT  is  satisfied  and  tnat  Eqs.  (3a)  tnrough  (3c)  are  equivalent.  We  caution, 
however,  that  the  proof  of  Eq.  (27)  in  Ref.  9  die  not  include  a  unitarity  cu¬ 
toff  for  strong  collisions.  One  must  be  careful  in  the  application  of  a  cu¬ 
toff  if  the  FDT  is  still  to  be  satisfied  (see  Section  5). 


By  using  symmetry  relations  (21),  (28),  anc  Eq.  (26)  in  the  form 

*fj(u  ♦  wf|)  -  +  -  a>fj)  (30) 

one  may  rewrite  Eq.  (18)  as  a  restricted  sum  over  pairs  of  states  with 
tf  i.e.  with  w^>0;  the  result  is 

<*Kw)  -  v  £  (31) 

iiif 

Cf>€j 


-  wfi>  *  '  wfi> 


The  above  result  has  the  nice  advantage  that  we  can  consider  a  total  contribu¬ 
tion  to  absorption  for  the  pair  of  state  t  . ,  e^.  We  shall  refer  to  the  first 
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term  in  Eq .  (32)  as  tne  posi ti ve  oscillator  term,  and  the  second  term  will  be 
reserved  to  as  the  negati  ve  oscillator  term,  i.e.  it  is  resonant  for'  w=-wf..<0. 
Note  tnat  computation  of  the  negative  oscillator  term  oniy  involves  tne  simole 
’■eolacement  of  arguments  Therefore,  for  the  most  part,  our 
discussion  will  concentrate  on  the  positive  oscillator  term. 


The  positive  oscillator  term  can  be  put  into  the  form  of  an  of f-resonance 
Lorentzian  with  a  frequency  dependent  halfwidth. 


where 
r  =N 


per 
ergs . 


n  =N 
rad 


to 


per 


rad 

i  s 


1  r<«-  uj^j) 

-  “«>  -  "rad  — 

/H  is  the  radiator  density,  and  where  the 
contained  in  the  halfwidth  function  r(u-w^i 


(33) 

perturber  aensity 
),  which  has  units 


We  now  write  down  the  equations  which  determine  the  halfwidth  function  r. 
^'•st,  we  do  this  in  a  simp!  1  fied,  schematic  fashion  which  suppresses  magnetic 
quantum  numbers,  Clebsch-Gordan  coefficients,  and  m-summations .  At  the  con¬ 
clusion  of  this  section,  we  shall  present  the  general  formulas  for  T  which  re- 
ouire  Clebsch-Gordan  reduction  (see  Appendix  C). 


The  simpl i fied  equation  for  r(w-w^.)  i: 


f(w  -  wfj)  a  2  e(Cj)e(€k)  2^  !  «X\V\0>  1 2 .  *5(€a0  +  hw  -  hwfj) 


W 

/3h(w-u>fi)  •  v-' 

♦  e  Nper  L  G(ej)G(€k) 


Jk 


•  2  |  <a'jVi/£?'>  i 

IfJ'k' 

where  for  simplicity. 


/?'>  iJ  T  5(€a,^,  -  hw  +  h«jj) 


(34) 


A 

JC 

n 

A 

a 

( 35a ) 

I0>  -  !iH'k'>  . 

(  25d ) 

|a'>  -  !jjJk>  , 

( 35c) 

!£'>  -  ijjJ'k'>  , 

( 35d) 

and  £a&=W  £a 1  S 1  =£a '  ”£  B '  " 

With  regard  to  Eg.  (34)  we  note  that  trie  contriduti ons  from  tne  initial 
and  final  states  (i.f)  of  tne  radiative  transition  ao  not  enter  in  a  symmetric 
fashion.  It  is  precisely  this  asymmetry  which  causes  tne  FDT  to  De  satisfied. 
At  line  center,  (w“Wr.)=0,  this  asymmetry  disappears.  Tne  nalfwidth  fjnction 
satisfies  the  same  symmetry  relation  as  ?  in  Eq .  (26),  namely  (at  the  risk  of 
introducing  more  indices),  if  we  define  r^(z)=r(z),  then 

ry-zi .  .  rlf(z)  <«> 

for  general  frequency  z  (which  is  to  be  substituted  for  (w-u>fi)  in  Eq .  (34)). 

The  first  (final  state)  contribution  to  the  halfwidth  function  has  a  sim¬ 
ple  interpretation ,  namely  it  is  essentially  Fermi's  golden  rule  averaged  over 
initial  states  and  summed  over  final  states.  The  interpretation  of  the  second 
term  in  Eq.  (34)  is  complicated  by  the  FDT  asymmetry.  Note  that  the  asymmetry 
involves  two  things: 

1.  the  temperature  dependenct  factor  exp[Sti(w-w^ . )]  in  the  second  term  of 
Eq.  (34);  and 

2.  the  off-resonance  energy  Ti(u-w^.)  enters  the  delta  functions  with  a 
different  sign  in  the  two  terms  of  Eq.  (34). 
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Genera i  Results:  We  close  this  section  by  writing  down  the  genera' 
eroression  for  the  halfwidth  ^unction,  including  C 1 eDsch-Gordan  coef*icierts 
ana  magnetic  auantum  numoe^s.  As  a  preliminary  to  listing  the  relevant  eaua- 
t;on,  we  note  frcm  £q.  (34)  that  we  wi1!  eventually  have  to  consider  a  summa¬ 
tion  over  t^ar'at'onal  states  of  the  form 

sum  *  Nn«r  2  e(ew)  <klV  k*>  <k'; V  k>  (37) 

p  kk' 

wne-'e  denotes  any  comolete  set  of  relative  translational  states. 


I”  Appendix  3  we  snow  that  this  sum  can  be  rewritten  as  (the  sum  is  sti 1 1 
an  operator  in  the  space  of  the  internal  molecular  coordinates) 


„ .  JWm*  £  £ 

2  '  T  '  (rrtf 

oo  OO 


dk'<R 


kf I v  I  Rk'f'mf ,  >  <  Rk'f ,  i v I  Rk'f'mf , : 


(38) 


In  tne  above  eauation: 

1  Rkfm/rl  ■  <r:Rkfmf> 

m  0.  <39> 

where  the  radial  wavefunction  (r  is  the  magnitude  of  the  relative  coor¬ 
dinate)  has  the  asymptotic  nonma’ i zati on 

r)  sin  (kr  -  Vi  It  +  fy(k)) 

2.  The  meaning  of  the  matrix  element  is 

<Rkfmf|V;Rk,£,mff> 

00  1  00 

*  1  <V  f  CJ(COS^)  j  r2drRkfm/r)  v(r>Rk'f'm*,<r)  (41) 

0-1  0  1 
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3.  Tne  quantity  z  is  a  translational  partition  sum  given  oy 


z  = 


-0«r 


~  { d*  e  ""k 

T\3  J 


-fee- 


(2t)3 


( 


m  \  3/2 


2x/3h2  / 


(42) 


ror  tne  nalfwidtn  function  we  obtain 


00 


r<“  -  «w  -  (-^)2  Se(ej)  J  - 

J  0 


-fc. 


00 


J  dk-E  £ 


{sum  op  ♦  sum  Gf)} 


(43) 


0  J'  If 

Here,  including  all  m-summations  which  must  De  reduced  using  Clebscn-Gordan 
algebra. 


sum  Of)  =  S  Tr^ff,  +  hu  _  hWfj) 
jf” 


E  E  E  E 


+  ^  m^m  mfmf,mf„  MM'  m^m^. 


Smj.mj,  Of  1mfmjjf  l^m^Of  1mf,m|jf  IjjiTij,) 


•  <jfmf,;  JM;  Rktfm£|V|jf«nif„;  J'M’;  Rk'£'m^,> 


.  <jf»mf„;  J'M';  Rk^m  i 


Ir  tne  celta  ^unction  expressions  above,  we  have  defined 

«ff"  *  («i,  -  %,)  *  ((J  -  fJ'>  +  <fk  -  <V>  •  C6) 

fii"  =  <€ii  -  %J  *  <«J  -  eJ>  *  <ek  *  {k'>  '  (47) 

^he  three  terms  on  the  right  hand  side  of  Eqs.  (46)  and  (47)  are,  respec¬ 

tively,  change  in  energy  of  the  radiator  (internal  energy),  the  perturber  (in¬ 
terna1  energy),  and  relative  translational  motion  in  a  collision  whe^e  f-*f"  or 

■;  -*■(  ii 

Some  details  of  the  Clebsch-Gordan  reduction  of  Eqs.  (44)  and  (45)  are 

absented  in  Appendix  C  eor  an  anisotropic  potential  corresoondi ng  to  a  di¬ 

pole-dipole  or  di pole-auadrupol e  interaction.  From  Appendix  C,  the  result  *or 
Eos.  (44)  and  (45)  becomes  (for  dipole-dipole  case) 
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sum  (jf)  *  w  x6<e^,  +  ho)  -  h<jjfj) 

V 

.  (2f  +  1)  (2J  +  1)  (f2  00s^2f'0)2 

w 

i<Rk^i  ;Rk'f'>!2 

•  I  <  jf  I  iMi  1  iif*  >  \2  1  <J!  !/*2 1  <J>  2  • 


(48) 
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At  a  convenient  point,  in  going  between  Eos.  (44)  and  (45)  and  Eqs.  (46)  and 
(49).  we  have  changed  names  of  dummy  variables  ,  ijH-j • i . 

Finally,  we  remark  that  the  two  terms  in  Eq.  (43)  correspond  to  Ander¬ 
son's  S^oute'"')  terms.  In  certain  problems,  one  has  to  also  worry  about  An- 
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derson's  $2(inner).  '  This  term  obtains  from  an  uncoupled  line  approxima¬ 
tion  to  Eo.  (51)  in  Ref.  9.  This  term  has  been  ignored  for  the  problem  of 
Drimary  interest  here  (H^O  se1 f-broadening)  because  it  involves  diagonal  ma¬ 
trix  elements  of  the  dipole  moment  operator,  for  the  case  of  dipole-dipole  in¬ 
teractions  . 

In  the  next  section  we  discuss  the  above  equations  further  in  terms  of 
quantities  manipulated  by  our  computer  codes.  In  doing  this,  we  also  identify 
the  ana’og  of  Anderson's  functions.  This  is  important  for  considerations 
o*  an  Anderson  cutoff  scheme  in  the  strong  collision  limit. 

3.0  REDUCTION  TO  COMPUTATIONAL  FORM 

The  Schroedinger  Eq .  (15)  which  determines  the  eigenstates  0 ^  ( k. ; r )  in 

Ecs.  (39),  (40),  (50),  and  (51)  has  been  solved  numerically  for  the  case  of  a 
Lennard-Jones  potential  with  wel1  minimum  energy  Es«  occuring  at  radius  r 
Some  details  of  the  matrix  element  calculations  will  be  presented  in  the  next 
section.  The  parameter  r  is  used  as  the  length  scale  unit  in  these  calcula¬ 
tions.  In  particular,  for  the  dipole-dipole  case,  matrix  elements  of  r  ^  are 


computed  as 


where  the  reduced  length  is  z-r/r  .  We  may  also  aenote  tne  aoove  matrix  eie- 

2  2 

ment  by  I(le;llc‘),  wnere  k.  /2m  's  tne  energy  asscciatec  witr, 
wavevector  k.  We  define  dioole  line  strengths  Dy 


l<j|iM!!i'>!2  =  M2Sr  .  f=* 

Here  u  is  the  scalar  dipole  moment  (u=l.S5xlC  esu-cm  is  tne  permanent  mo¬ 
ment  for  H^O).  and  the  line  strengths  satisfy  the  sum  »*ole 

E  S:.,  *  1  .  (54) 

r  1 

For  programming  convenience  we  will  eventually  (but  not  immeciately)  co"- 

_  i 

vert  all  energies  from  ergs  to  cm  *  using 


€(ergs)  »  hce(cnr’)  . 


For  the  case  of  sel f-oroadeni ng  through  dipole-dipole  interactions,  we  compute 
and  store  the  following  dimensionless  (squared)  matrix  element 


Mf(c;  €'>  -  (~— ~)2  C(€)c(0 

S 

•  E  (f200|f2f'0)*!l(fe;  f'e*)la 


In  tnis  expression,  with  e  in  cm 


I  8mcr|  \  V2 
V  he  ' 


is  a  density  of  states  factor  with  units=cm.  That  the  M-matrix  element  aoove 
is  dimensionless  follows  from  the  fact  that 


/  2  \ 2 

te?) 


has  units  cm  , 


while  p(e)p(e')  has  units  cm4,  and  everything  else  is  dimensionless. 
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We  next  introduce  a  function  which  plays  a  role  analogous  to  Andersen's 
1 2  13 

s^  function.*  ’  We  write  tne  function  ,-n  the  *orm 

S20t(£,  J,  t;  o)  -  ojfj)  =  sj  (e,  J,  f;  w  -  o;fj) 


0h(<i)-a)f)  j 

+  e  S2(€,  J,  f;  o)  -  ojfj) 

(53) 

f  * 

For  the  dioole-dipc'e  case  the  (d'mens'orless)  functions  s,  ana  s, 

are  g’ven 

py 

S2  (€,  J,  f;  (j)  -  oofj) 

=  E  Sjfjf',  Sj,,  M^[e;  e  +  e\Qt  +  h(w  -  wfj)J  , 

J  jf.  J' 

(59) 

S2  (€,  J,  f;  oj  -  ofj) 

=  J)  Sjjj,',  Sjj,  M<ce  e  +  4oi  -  h<“  -  "(ill  • 

V  J' 

(60) 

with  fr0(  *  (fjf  ~  ^jf')  +  (Cj  ~  ej')  • 

(61) 

erot  -  +  <€J  "  €J’)  • 

(62) 

That  the  aDove  plays  the  role  of  Anderson's  s^^function  (at  least  in 
tne  line  center  limit  w-w^.-'O)  will  presently  be  demonstrated.  The  FDT  asym¬ 
metry  detween  the  final  and  initial  state  contributions  is  again  evident  in 
Eos.  (58)  through  (60). 


For  the  halfw’dth  'unction  of  Eg.  (45),  we  may  now  write 


r<w-  wfj) « 


L 


J 


e(€j)  (2J  +  1)  a  (e,  J) . 


(63) 
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In  tne  aoove  exoression.  o(e.J)  is  a  cross  section  ( cm“ )  g’ven  oy 

o(€.  J)  -  (-|~-)  E  (2f  +  1)  Sa0t(e.  J.  *;  «  -  wfi)  •  ( 6i ) 

Before  going  on,  we  note,  if  Eqs.  (63)  and  (64-)  are  comDinea.  tnat  the  rasters 

7  ... 

(irfi  /2me)  cancel.  ihis  is  userul  in  terms  or  writing  tne  results  in  alterna¬ 
tive  forms. 

Next,  in  tne  limit  of  large  1,  ana  using  l=kD.  where  b  is  the  class ical 

2  2 

impact  parameter,  with  e=fi  k  /2m,  one  can  see  that  o  in  Eq.  (70)  takes  the  ap¬ 
proximate  form 

00 

o{€,  J)  s  |  2t  bdb  s|ot(6,  J,  b;  w  -  o>fj)  ( 65 ) 

This  is  the  usual  form  which  relates  the  s^  function  to  tne  cross  section  in 
Anderson's  formalism. 


To  complete  the  analogy  with  Anderson's  theory  we  rewrite  the  average 
over  energy  in  Eq .  (63)  in  terms  of  a  Maxwel 1-Boltzmann  velocity  average.  We 
use 


(66) 
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r  (w  -  wfj)  =  r(«  -  ojfi) 


(68) 


7  is  the  halfwidth  furction  in  '■ad/sec.  we  obtain 

00 

r[w  -  o>fj)  =  npef  j  Vdvf(v) 

0 

•  E  e(«j)(2J  +  1)  a(v,  J) 

with  J 


(69) 


f(v)  = 


v2  e 


-0  J  mv2 


00 


I 


v2  e 


-0  -J-  mv2 


dv 


(70) 


tre  Maxwe  ’  '  -3c '  tcm.ann  ci  sfi  but;  on  function.  We  can  rewrite  this 

F(0J  -  Wfi)  =  nper  E  e(€j)(2J  +  1)  <v  o(v.  J)>  .  (71) 

J 
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>  the  line- center  '  ’  m  i  t  this  has  the  form  of  an  Anderson  ’  type  theory, 
'rc'ucmg  velocity  averaging.  With  an  aporopriate  definition  of  the  average 
c^oss  section,  Eo .  (71)  can  be  exDressed  as 


r  .  nper  v(T)  a(T)  . 


(72) 


wne^e  v  is  tne  mean  thermal  velocity. 


"h's  essentially  completes  the  analogy  of  the  present  theory  with  Ander¬ 
son's  forma 1 sm .  Certain  questions  remain  as  to  wnetner  Anderson's  cuto** 
procedures  for  s^  can  be  tax.en  directly  over  for  the  *ar-wing  problem. 

finally,  we  wish  to  write  down  the  previous  equations  in  the  form  in 
which  tney  were  actually  programmed  for  calculation.  This  basically  involves 
combining  Eqs.  (63)  and  (64),  canceiling  the  factors  of  (irfi<V2me ) ,  and  con¬ 
verting  all  energies  to  cm  including  the  partition  sum  denominator  inte¬ 
gral,  which  can  be  written 
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_  3  _ 

X  e 


f  at- 


nere  e  .  e  are  in  cm  anc  we  nave  definec  some  average  tnermal  Quantities  a: 
fo:lows:  the  Dasic  definition  is  tne  mean  tnermal  velocity 


(+-¥) 


In  terms  of  v  we  also  define 


€  (ergs) 


2  m  v 2  »  hce(cm'’)  , 


i.e.,  €(cm*’) 


(a  kBT\ 

It  he  J 


We  caution  tnat  Eg.  (75)  is  simply  a  def i ni tion  of  t,  and  I  is  not  the  average 
energy  in  a  thermoaynami c  sense.  Finally,  we  define  an  average  DeBroglie 
wavelength  by 


JL  ,  h  \  Va 

mv  V  2mcT(cnv1)  / 

Equations  (33),  (63),  and  (64)  can  now  be  re-expressed  as 

♦if"  -  “fi>  -  nr,d  —  _L 

rad  he  T  (W  -  Wfj)2 


n«  -  «fj) 


X’T  (  dte-'3h« 

0_ _ 

00  ~ 

J  de'  /-4r-  e~^hc€ 


E  G(Cj)(2J  +  1)  £  (2f  ♦  1)  s|0t(€,  J.  (■  CJ  - 
J  ( 
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In  the  above  equations,  all  e's,  w's,  and  T  are  now  in  cm  The  units  of  r 
are  manifest,  namely  everything  is  dimensionless,  except 

^nper  X3)  e  -  €  »  cm*1  . 

To  conclude  this  section,  since  it  would  be  totally  prohibitive  to  actu¬ 
ally  calculate  the  absorption  coefficient  (by  summing  the  contribution  from 
all  Mnes),  we  have  to  proceed  in  a  different  fashion.  For  certain  experimen¬ 
tal  cases  it  is  known  what  the  simple  impact  approximation  predicts  for  ab¬ 
sorption.  This  ’  s  subject  to  the  difficulty  that  the  impact  approximation 
does  not  satisfy  the  FDT,  so  that  the  impact  results  depend  on  which  formula 
(Eqs.  (3a)  through  (3c))  one  starts  from. 

We  shall  make  comparison  with  the  so-called  "hyperbolic  tangent"  impact 
formula.  This  formula  has  some  advantages,  and  certain  predictions  from  the 
formula  are  available. 

Before  writing  out  tne  impact  formula,  we  combine  Eqs.  (3b),  (31),  (32), 
and  (78)  to  deduce  that  the  contribution  to  x"  associated  with  the  pair  of  en¬ 
ergy  levels  is  proportional  to 


X"<«»FOT  *  0  -  e-0hc“) 


r(M  ~  “V  +  r(~M  ~  > 

(W  -  Wfj)2  (-0)  -  OJfi)2 


.  (80) 


wnere  the  two  terms  in  the  bracket  are  the  positive  and  negative  oscillator 


terms,  respectively. 


The  corresponding  "hyperbolic  tangent"  impact  formula  is 
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j_(0>  -  Uff  (-0J  -  u)fj)2  J  • 

In  writing  down  the  above  results  we  have  dropped  some  constant  factors  wnich 
will  cancel  in  the  ratio 

R(w)  »  X"Mfdt  /  X"(co)jmpact  .  (32) 

One  of  tne  constant  factors  which  has  been  dropped  is  the  radiator  cen- 

sity,  n  As  a  matter  of  fact,  the  perturber  density,  n  ,  also  cancels  ’  r 

rad  per 

the  ratio  Decause  alj  of  the  halfwidth  functions  in  Eqs.  (80)  and  (81)  are 
linear  in  n  .  Thus  R(w)  is  independent  of  both  radiator  ana  perturber  aen- 
sity. 

The  only  remaining  density  dependence  of  interest  is  in  the  halfwidth 
function  itself.  We  have  done  all  calculations  at  a  fixed  density  correspond¬ 
ing  to  the  density  at  one  atmospnere  pressure  and  room  temperature,  i.e.  we 
ta*e  for  the  perturber  density  the  fixed  value 

np^  •  V273/296)  •  (83) 

with  n^toschmidts  number.  R.H.  Tipp;rg  nas  referred  to  the  quantity  above  as 
a  "hot  Loschmidt" . 

To  concluae  this  section,  we  note  as  ,  wnere  the  negative  oscillator 

terms  in  Eqs.  (80)  ana  (81)  become  negligible,  that  the  two  formulas  approach 
each  other.  Similarly,  the  formulas  coalesce  at  the  negative  oscillator  reso¬ 
nance  Of  course,  both  formulas  are  divergent  at  the  line  center. 


4.0  CALCULATION  OF  MATRIX  ELEMENTS 


I n  terms  of  computation  time,  calculation  of  the  M-matrix  elements  in  Eq . 
(56)  is  oasica'iy  governed  by  the  speed  with  which  the  simple  I-matrix  ele- 
merts  in  Eq.  (52)  can  be  generated.  We  nave  attempted  to  mane  this  process  as 
efficient  as  possible. 

To  discuss  simp! '  fi cat’ ons ,  we  begin  by  rewriting  Eq .  (52)  as 

00 

\(le.  (’€•)•  J  d Z0f(e;Z)  <t>( Z) ,  (84) 

where  e.  z 1  are  the  energies.  In  considering  the  above  integral  we  remark, 
fo,‘  most  cases  of  interest,  that  1,1'  and  e,c'  are  "not  all  that  different". 

To  be  more  precise:  for  the  dipole-dipole  case,  the  Clebsch-Gordan  coef- 
ci ci ent  (  1200  121 1 0 )  in  Eq .  (56)  leads  to  the  selection  rules  ( 1 ' - 1 )=C . ±2 . 

rhe  statement  in  regard  to  e.e'  is  simply  that  if  e-e1  is  large,  the  matrix 

elements  tend  to  be  small.  The  point  of  the  above  comment  is  simply  that,  for 

most  cases  of  interest,  the  wavefunctions  in  the  integrand  are  such  that  they 

can  generally  be  calculated  using  the  same  approximations. 

With  the  aDove  remark  in  mind,  there  are  essentially  three  situations  of 
interest  which  depend  on  1,  e,  and  the  total  isotropic  potential 

W  f)  -  vu<r>  ♦  (85> 

In  the  above  expression  the  first  term  is  the  Lennard-Jones  potential  and  the 
second  term  is  the  repulsive  centripetal  potential,  which  is  clearly  '-depen¬ 
dent. 


In  Fig.  1  we  show  a  schematic  Vg(r,l)  potential,  along  with  three  energy 
levels  e,,  e,,  and  e,.  We  will  provide  a  few  details  at  the  end  of  this  sec- 


tion.  but  briefly  the  wavefunctions  associated  with  tne  above  energies  a^e 
calculated  as  follows: 

Case  ( 1 ) .  Energy  level  (sufficiently  low  e.  sufficiently  r,:g-.  1;.  we 
call  this  tne  centri petal  regime.  Classically,  tne  turning  pcirt  is  sucn  tnat 
the  wavefunction  never  gets  in  to  "see"  the  Lennard-Jones  potential.  Tne  wa¬ 
vefunctions  are  then  approximately  eigenstates  of  a  Hamiltonian  witn  .=0. 
The  approximate  eigenstates  are  Bessel  functions  anc  tne  integral  in  Eq.  (34) 
can  be  expressed  in  terms  of  the  hypergeometri c  function. 

Case  (2).  Energy  level  z^  ■  We  call  this  tne  intermeal  ate  region.  In 
this  case  the  wavefunctions  <t>  are  obtained  numerically,  as  is  tne  integral 
(84). 

Case  (3).  Energy  level  (sufficiently  high  e,  low  1).  We  cal:  this 

the  hard  core  regime.  Here  the  wavefunctions  relatively  quickly  assume  tneir 

asymptotic  form  (Eq.  (40))  close  to  the  turning  point.  The  procedure  in  this 

case  is  to  use  numerical  solutions  out  to  a  radius  Z  =R  /r  where  it  nas 

max  max  s 

been  determined  that  both  $’s  have  assumed  (to  a  good  approximation)  their  as¬ 
ymptotic  forms.  Beyond  this  radius,  tne  asymptotic  form  (40)  is  employed,  ana 
tne  resulting  integral  can  be  integrated  in  terms  of  sine  and  cosine  inte~ 
gral s . 

We  close  this  section  with  a  few  details  of  the  calculations. 


Case  ( 1 ) .  Centripetal  Regime.  The  approximate  wavefunctions  are 


0f(«:  Z)  -  (KZ)  if  (KZ)  , 


(86) 


V2  (KZ)’*  Jf  +  Vi(KZ). 


(87) 
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re- 


where  K=k.rs>  and  where  j,  J  a<-e  spherical  and  ordinary  Bessel  functions, 
spectively.  The  integral  (52)  becomes 

Kff;  Ce’)  *  -|-  (KK')Vs 


T  h  i  s 
page  497, 


oo 

•  J  dZ  Je  +  v,!*2*  Jf'+ %<K,Z>  <  as  > 

is  a  Weber-Schafhei tl  in  type  integral  (see  Abramowitz  and  Stegun^, 
Eq.  11.4.34).  The  result  can  pe  written 


(KK')V2 


K 


t  +  Va 


r«  +  -y-  ) 


(K7 


£-Va 


r<|  +  +  {> 


where  Al=(l'-1)=0,  ±2,  T  is  the  Gamma  function,  and  where  F  is  the  Gauss  hy¬ 
pergeometric  functior.  The  above  formula  is  valid  only  for  K/K1  <1;  however, 
Eq.  (88)  is  seen  to  be  invariant  under  K i  K'  and  1  £  1 '  ,  so  matrix  elements 
w’th  K/K 1 > 1  can  easily  be  obtained  from  symmetry.  We  note  that  the  function 
in  Eq.  (89)  has  the  form  F(a,b;c;z)  with  c=a-*-b-*rn,  and  with  m=2.  For  this  ar¬ 
rangement  of  arguments  the  F-function  is  most  conveniently  calculated  from  Eq . 

14 

15.3.11,  page  559  of  Abramowitz  and  Stegun. 

Case  (2).  Intermediate  Region,  The  0's  are  solved  for  numerically  via 
Numerov's  method  using  a  computer  code  developed  by  R.J.  LeRoy.^  The  integral 
(84)  is  evaluated  numerically. 


Case  ( 3 ) ,  Hard  Core  Regime.  We  use  tne  same  procedure  as  in  Case  (2) 


aDove  out  to  a  radius  Z  ,  then 

max 


max 

J  -I 


The  second  integral  above,  for  the  dipole-dipole  case,  is  of  tne  form 


f  dZ 

I2  -  J  —  sin(KZ  +  A)  sin(K'Z  ♦  B) 


wnere  the  pnases  A.B  nave  been  determined  from  the  Z<Z  region.  By  repeated 

rnd  x 

integration  by  parts,  Eq.  (91)  can  be  expressed  in  terms  of  sine-cosine  ;nte- 
grals.  We  define  K±=K±K‘  and  C±=A±B.  Then 


l2  -  I  -  I 


« - r  cos<K-2max  *  C-> 

-  I T~—  sin(K-  Zmax  +  c-> 


♦  —  { cos(C  J  Ci  ( !  K_  |  Zmax) 

-  sign  (KJ  sin  (CJ  si  ( |  ,  Zma^) } 


here  sign  (K_)=tl  depending  on  whether  K_  is  positive  or  negative.  In  Eq. 
(93),  si  and  Ci  are  the  sine,  cosine  integrals  given  by  (see  Ref.  14,  pages 


231-232) 


-j-f 


sin(t)  , 


-J-f 


cos(t) . 
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The  integral  I+  in  Eq .  (92)  is  obtained  by  replacing  C_-*C+  and  K_-*K+. 

The  I  and  M-matrix  elements  were  computed  using  the  aoove  methods,  and 
the  M-matrix  elements  were  stored  in  a  loox-up  table.  Wnen  M-matrix  elements 
are  required  (Eqs.  (59)  and  (60))  at  nonstored  e.e'  values  (as  is  virtually 
always  the  case),  interpolation  from  the  look-up  table  is  employed. 


5.0  RESULTS  AND  DISCUSSION 

In  the  previous  sections  we  have  presented  a  far-wing  lineshape  theory, 
assuming  dipole-dipole  collisions  to  be  the  dominant  broadening  mechanism. 
For  close  collisions  one  might  expect  the  higher-order  multipole  interactions 
to  become  important.  The  selection  rules.  AJ=0,±1,±2.  for  the  dipole-quadru- 
pole  interaction  allow  for  larger  changes  in  rotational  energy  in  Eqs.  (46) 
and  (47)  than  oDtained  from  the  dipole-dipole  selection  rules.  AJ=0,±1.  Tlv's. 

; n  turn,  makes  it  easier  to  satisfy  the  delta  functions  for  conservation  of 
energy  in  Eq$.  (34),  (44),  and  (45)  when  (w-w^.)  Decomes  large. 

We  have,  in  our  calculations,  explicitly  included  the  di pol e-quadrupol e 
and  quad^upole-dipole  interactions.  (The  distinction  being  radi ator-perturber 
and  perturbe'— radiator,  respectively.)  In  the  case  of  the  quadrupole-dipole 
interaction  there  is  also  a  contribution  associated  with  Anderson's 
S7(  i  n-'e’-) .  ’ 1  This  contribution  has  been  ignored  in  the  present  calcula¬ 
tions. 

The  essential  modifications  necessary  to  treat  the  dipole-quadrupole  and 
quadruool e-di pol e  interactions  are  given  in  Appendix  C.  Somewhat  more  explic- 
it’y,  tne  modifications  can  be  simply  expressed  as  presented  below. 
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For  the  di pol e-quadrupol e  case,  for  the  functions  in  Eos.  (59)  anc 
(60),  the  following  substitution  is  to  be  made: 

2 

4tt  1 S  2  7 

—  Si  i  SJJ'  ^  (T)  (2?r  )Si  i  St-*  (96' 

In  this  suosti tut i on ,  tne  difference  in  numerical  factors  arises  precise', 
from  Eq.  (C30a).  Also  in  tne  above  equation.  0.,,  is  a  ouacrupc’e  strenotr, 

UVJ 

defined  in  terms  of  tne  reauced  matrix  element  (wnicn  appears  ;n  Ecs.  (C27) 
and  ( C 29 ) )  Dy 

Z:<JiiQ!iJ’>!2  *  q2  £  St  “  q2-  (57) 

j.  J(  jj 

2 

That  is,  q  is  tne  square  of  the  scalar  quadrupole  moment  (it  is  ircepencent 
of  J),  and  tne  quadrupole  strengths  satisfy  the  sum  role 

(98) 

For  the  quadrupole-dipole  case,  the  substitution  analogous  to  Eq .  (96)  is 


*r  sj  ,  S,j,-><£)2<2*2)  G.  .  s„,  (99) 

JfJf.  JJ  7  JfJf,  JJ*. 

The  dipole  strengths  Sjj,  were  computed  using  the  asymmetric  rotor  codes  de¬ 
veloped  at  AFGL.  The  quadrupole  strengths  G,,,  were  computed  using  a  program 

<JvJ 

developed  by  one  of  the  present  Authors,^  which  can  De  attacned  as  an  upcate 
to  the  AFGL  package. 


The  remaining  modifications  are  to  the  translational  matrix  elements  and 
tnese  are  identical  for  the  di pol e-quadrupole  and  quadrupole-dipole  cases: 

first,  the  factor  z  ^  in  the  I-matrix  element  (Eqs.  (52)  and  (34))  gets  re- 

-4 

placed  by  z  ;  then  in  Eq.  (56)  for  the  M-matrix  element,  the  following  sub¬ 
stitution  is  necessary: 
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zra'  (Eas.  (9C)  through  (35) ) .ca'culations  discussed  i "  Sect'or  4  tnerp  ' s  Bi 

-3  -4  .  .  . 

cbv’cus  mcc'  *  ■  cat'  or  cue  to  z  -z  ,r  tne  .-matrix  element. 


We  ncw  oucte  tne  soec ’  *  i  c  parameters  upon  which  ow  presert  za 1 cul ations 
are  cased.  -or  toe  .er-arp-Jones  coterfa’  we  nave  cnoser  the  parameter  va'~ 
^es,*  E  =243  cm  '=potent':a'  nrr'mum,  anc  '•_  =  2.C59  cm=radius  at  wel’  nin’mum. 

3  3 

the  oer'nanert  c;Do’e  moment  we  rave  usee  the  we" '-known  value 

-  *  Q 

u=l.S5xlC  esj-cm=i.85  Debye.  ►o’’  tne  H^C  Denmanent  ouadnuDoie  moment  we 

18  .  . 

-ave  emc  loved  the  review  a-tic'e  o*  ^trocryn  and  Strogryn  to  inte^  an  aver- 

2s  2 

age  measured  va'ue  cr  o=1.158xiC  “  esu-cm  . 

The  cn^  remaining  adjustable  "parameter"  in  the  tneory  is  the  manner  in 
w-'-’ch  a  'jn'ta-ity  cutof*  ;  s  emp'oyed  in  tne  strong  collision  (small  imoact  pa¬ 
rameters,  classically)  limit.  cor  sufficiently  small  (w-w^.),  there  1-s  no 
cc-bt  that  such  a  cutoff  is  necessary,  because,  the  cutoff  is  eliminated, 
the  resulting  halfwidtrs  at  the  line  center  are  large  by  a  factor  of  order  40 
comoared  to  experiment.  Sufficiently  far  in  the  wings  the  cutoff  can  be  drop¬ 
ped.  However,  exactly  how  one  should  interpolate  between  these  two  'Suits  re¬ 
mains  somewhat  of  an  open  Question. 

4noersor's  simplest  cutof*  procedure  as  aool’ed  to  ’is  to  reset 

STct=l,  '  e  S~ot  (as  calculated  by  Eq.  (58))  is  such  that  it  exceeds  j"ity. 

C.  C. 

An  alternative  (and  somewhat  smoother)  exponentia'  form  of  the  cuto£f  is  to 
replace  S0  by 

(S20t)eff  "  1  ~  e_S2  •  (101) 


31 


Theoretical  justification  for  this  type  of  cutoff  may  be  found  in  tne 
literature"9'22.  From  Eq .  (101).  if  S^^l.  we  obtain  ( =  5^ . 

I",  on  the  other  hand.  S^ot»l.  Eq .  (101)  reauces  to  ( S*01)^  c  =  i . 

The  above  cutoff  schemes,  based  on  S^0^,  suffer  from  a  major  problem  in 
the  far-wing  case.  Due  to  the  exponential  factor  in  tne  second  term  of  Eq . 
(58),  the  FDT  cannot  be  satisfied  unless  the  cutoff  is  applied  separately  tc 
s!  and  s\  in  Eq .  (58). 


In  the  absence  of  a  rigorous  prescri pti on  for  cutoff,  we  have  emoloyec 
tne  following  methods  which  satisfy  FTD  anc  yield  an  which  is  bounded  by 
unity  as  in  Anderson's  theory.  The  first  metnod  is  simply  to  reset  (cr 
S2)  equal  to  i/2  if  S2  (or  S^)  exceeds  1/2.  An  exponential  form  of  this 


cutoff  scheme  is 


(102) 


Again  for  s£»l  this  yields  (s|)eff=sj,  while  for  s£«l.  we  obtain 
1/2. 


In  Tables  1  through  4  we  present  tne  results  of  our  calculations  of  I^C 
far-wing  absorption.  The  calculations  have  been  carried  out  using  the  expo¬ 
nential  form  (£q.  (102))  for  the  separate  cutoffs  on  and  S^. 

Tne  calculations  have  been  carried  out  for  tne  strong  pure  rotational 
transition  4, 1,4-5, 0,5  with  line  center  frequency  at  wfi=100.51  cm-1.  In  Ta¬ 
bles  1  and  2  we  present  results  for  a  photon  energy  w=500  cm"1.  This  means 
that  the  positive  and  negative  oscillators  are  "off  resonance"  by  approxi¬ 
mately  400  cm  *  and  600  cm  1 ,  respectively.  The  columns  in  the  tables  are  as 
fol lows . 
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Column  1  =  temperature  °K. 

~  -t 

Column  2  =  T(0),  the  line  cente*-  halfwidth  in  cm  A  as  given  Dy  Eq .  (79)  with 
( w-wf . )  set  to  zero . 

■<*w 

Column  3  =  r(w-wp. )=the  positive  oscillator  halfwidth  as  computed  from  Eq . 
(79). 

Column  i  =  :(-w-w*,)=negative  oscillator  halfwidth  computed  from  Eq .  (79). 

Column  5  =  x"(w)cn-  as  computed  from  Eq .  (80). 

Column  6  =  x"(u>)-  4  as  computed  from  Eq.  (81). 

x  mp  3i  c  v 

Column  7  =  x"(u)pg,/x"(<»>)  jmpact  as  computed  from  Eq .  (82). 

TaDle  1  presents  results  for  w=500  cm  *  in  which  only  dipole-dipole  con¬ 
tributors  are  included.  Table  2  is  for  w=500  cm  but  including  dipole-di¬ 
pole.  di pol e-auadrupol e .  and  ouadruoole-dipol e  interactions.  From  the  line 
center  halfwidths  in  Tables  1  and  2  we  note  that  the  addition  of  the  dipole- 
quadrupole  interaction  represents  about  a  2%  effect,  so  that  the  dipole-dipole 
interaction  is  strongly  dominant.  Fom  fixed  density,  the  increase  in  the 
halfwidth  with  temperature  is  mainly  associated  with  the  increase  in  v  which 
appears  in  Ea .  (72).  (See  Fig.  2  for  a  plot  of  the  temperature  dependence.) 
The  values  for  the  halfwidths  are  approximately  what  one  would  expect  from  Ar- 
derson's  theory  ^or  self  broadening  of  Iow-J  H^O  transitions. 

The  positive  and  negative  oscillator  halfwidths  are  seen  to  be  reduced 
relative  to  the  line  center  halfwidth  (the  impact  halfwidth)  due  to  the  diff’- 
culty  of  conserving  energy  far  off  resonance.  This  also  accounts  for  the  fact 


33 


that  the  negative  oscillator  halfwidtn  is  small  compared  to  the  positive  os- 
ci 1 lator  hal fwidth. 

It  is  very  interesting  to  note  that  the  addition  of  tne  dipole-quadrupcle 
and  quadrupole-dioole  interactions  in  Table  2  increases  the  positive  oscilla¬ 
tor  halfwidths  by  a  factor  of  order  2  and  increases  tne  negative  oscillator 
halfwidths  by  an  order  of  magnitude!  This  empnasizes  the  importance  of  tne 
nigner-order  multipole  interactions  for  the  far-wing  oroblem.  The  origin  of 
this  can  be  ratner  simply  stated:  the  translational  matrix  elements  ( Ea . 
(56))  tend  to  be  small  unless  tne  difference  in  translational  energies. 
t-z'  ,  -is  small,  typically  less  tnan  20  cm  V  This  being  the  case,  wren 

fi(u)-w^;)  is  large,  tne  only  thing  left  to  conserve  energy  is  the  cnange  in  ro¬ 
tational  energy.  Because  the  selection  rules  on  AJ  for  the  dipole-quadruple 

* 

and  higher-order  multipole  interactions  allow  for  progressively  larger  cnanges 
in  rotational  energy,  tnese  collisions  are  more  effective  in  broadening  tne 
wings. 

For  some  fixed  value  of  fi(  ww^. ) ,  there  may  De  some  particular  order  of 
multipole  interaction  which  dominates;  however,  tne  convergence  of  tne  multi¬ 
pole  series  cannot  be  assured.  Since  the  higher-order  multipole  interactions 
are  progressively  more  effective  at  close  separations,  this  also  empnasizes 
the  point  that  the  far  wings  are  dominated  by  strong,  short-ranged  collisions. 

The  results  for  the  ratio  R(w)  (Eq.  (82))  in  Column  7  indicate  too  low  a 
level  of  absorption.  At  1000  cm  it  is  known  that  the  impact  formula,  Eq. 
(81),  yields  an  absorption  which  is  high  compared  to  experiment  by  a  factor  of 
approximately  two.  Thus  at  <*>=  1000  cm  ",  one  would  1  i  x.e  to  find  the  ratio  R(w) 
=  1/2.  However,  at  w=500  cm  j,  R(w)  is  already  small  of  oraer,  R(w)=.20. 
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The  results  for  w=1000  cm  ^  are  presented  in  Tables  3  and  4.  These 
tables  show  the  same  general  trends  as  discussed  for  Tables  1  and  2.  The  ef¬ 
fect  of  adding  the  dipole-ouadrupole  interactions  is  even  more  pronouncec.  the 
>-atio  R(w)  is  raised  by  two  orders  of  magnitude  upon  inclusion  of  the  diooie- 
quadrupole  interactions.  In  Tables  1  and  4,  Column  7  shows  little  indication 
of  the  negative  temDerature  dependence  observed  experimentally. 

Future  Work. 

Altnougn  our  present  calculations  indicate  a  level  of  far-wing  absorption 
wh':ch  is  too  small  to  account  for  experimental  observations,  these  calcula¬ 
tions  nave  shed  light  on  a  numoer  of  imoortant  issues  which  remain  to  be  ad¬ 
dressed  . 

1.  A  more  rigorous  aDoroach  to  the  far-wing  unitarity  cutoff  procedure 
should  be  sought.  A  more  relaxed  cutoff  procedure  would  certainly  en¬ 
hance  the  far-wing  absorption. 

2.  The  sensitivity  of  tne  far-wing  absorption  to  the  lennard-Jones  parame¬ 
ters  should  be  investigated.  In  particular,  a  deeper  LJ  potential 
should  enhance  the  translational  matrix  elements. 

3.  It  would  be  desirable  to  extend  the  present  calculation  to  include 
auadruDole-quadrupol  e  collisions.  This  should  give  some  indication  of 
the  convergence  of  the  multipole  expansion.  If  divergence  is  indi¬ 
cated,  some  more  realistic  formulation  of  the  short-ranged  anisotropic 
interaction  will  be  required. 

4.  It  would  be  extremely  interesting  to  apply  the  present  formalism  to  the 
problem  of  broadening  of  H^O  in  the  far  wings.  It  is  known,  experi¬ 
mentally,  that  broadening  leads  to  a  greatly  reduced  continuum  ab- 
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sorption.  From  our  present  point  of  view  it  is  tempting  to  speculate 
that  this  reduction  is  related  to  the  small  rotational  constant 
(2  cm  for  which  in  turn  implies  much  smaller  changes  in  rota¬ 
tional  energy  in  the  collision  processes. 

Deep  resonant  and  true  bound  states  may  be  important,  including  tran¬ 
sitions  from  bound  to  continuum  states.  The  present  formulation  is  not 
applicable  to  the  bound  state  problem;  however,  the  singie-perturber 
approximation  whicn  concentrates  on  tne  interaction  between  a  single 
radiator  ana  perturoer  should  remain  valid. 
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LIST  OF  TABLES  AND  FIGURES 


TABLE  1:  Summary  of  absorption  for  w=500  cm  *  for  various  temperatures. 

Only  dipole-dipole  collisions  are  included. 

TABLE  2:  Summary  of  absorption  for  w=500  cm  ^  for  various  temperatures.  All 
three  interactions  (dipole-dipole,  di pole-quadrupole ,  quadrupole- 
dipole)  are  included. 

TABLE  3:  Summary  of  absorDtion  for  w=1000  cm  ^  for  various  temperatures. 

Only  dipole-dipole  collisions  are  included. 

TABLE  4:  Summary  of  absorption  for  w=1000  cm  ^  for  various  temperatures. 

All  three  interactions  (dipole-dipole,  dipole-quadrupole,  quadru- 
Dole-diDole)  are  included. 

FIGURE  1:  A  schematic  of  the  total  isotropic  potential,  Vg(r,1),  for  1>0. 

Three  values  of  the  translational  energy  e  are  indicated. 

FIGURE  2:  Temperature  dependence  of  the  linecenter  (impact)  halfwidths.  It 
is  assumed  that  a  power  low  of  the  form 

r(T)/r(296)  =  (T/296)m 

is  approximately  valid  so  that  the  quantities  plotted  are 
In  [r(T)/r(296)]  vs  In  [T/296 ]  . 

The  slope,  m,  from  the  figure  is  found  to  be  m  =  .365. 
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APPENDIX  A 


Derivation  of  The  Single-Perturber  Approximation  from  the 
Autocorrelation  Function  in  the  Time  Domain 


*r  th's  Appendix  we  derive  tne  si ngl e_perturDer  approx  *  mat-  or  the 

t 1 me~doma 1  r  treatment  o*  the  dipo’e  autocorrelation  funct'icr',  as  presentee  in 
Ref.  8.  crom  Ref.  3  we  make  use  of  Eqs.  (3.20b).  (3.25).  (4.13)  t-'rpugr 
(4.15),  and  (4.27),  (4.28).  ’’hese  equations  ignore  y-sums  and  ClePsco-Gordan 
coeff i ci ents .  and,  for  simplicity,  we  continue  this  omission. 


In  seconc  orcer  perturbation  theory  in  the  anisotropic  interaction,  the 
d;poie  autocorrelation  function  can  be  written 

"  ®*P  (NpeW)  +  W  +  Ki(d)]|  .  (Al) 

The  terms  Kr(t),  K,(t)  correspond  to  Anderson's  stouter)  terms,  while  K.^(t) 
corresponds  to  skinner).  We  shall  ignore  K.f(t)  for  the  reasons  given  at  the 
end  of  Section  2. 

~he  expressions  for  Kf,  K..  are 

W  *  -  llr  ^  C(€j)  C(ek)  £  |  < or | V ; 0 >  1 2  f(t)  ,  (A2) 

Jk  jf«J'k' 

W  "  "  TT  ?  G(«j)  C(€k)  2  <a'‘Vij3'>  2  f(t) .  (A3) 

J*  j-jv 

wnere,  for  simplicity, 


B>  = 


jfJk> 

,  J '  k’  > 


B’>  = 


j,Jk> 


j  7- ,  J  ‘  k  •  > 


We  also  define  fiu  =e  -e.  and  'flu  Iqi-e  i~eci. 

ab  a  3  a  S  a  b 

The  time  dependence  is  contained  in  tne  functions 

«*>  -  11  -  cos<“O0'»i 

-  -  s'n<"a/3,)l  ■  (Ai) 

'<•>  *  -^7  I1  -  e  0 

*  <-C/3'  Ia:“'3'('  "  l(3h)  "  e^0'*3'  sin(wa,^,t)l .  (AS) 

Tne  FDT  asymmetry  between  tne  final  Eq .  (A4)  and  initial  Eq .  (A5)  is  cor.ta^rec 
in  the  temperature  dependent  fjfiu  (Q,  factors  in  Eq .  ( A5 ) . 

a  p 

The  single  perturDen  approximation  assumes  tlYat  the  far-wings  will  De 
governed  by  times  sufficiently  small,  and  densities  sufficiently  low  trat  Eo . 
(Al)  may  be  approximated  by 

s  1  +  NpKf(t)  +  NpKj(t)  ( A6 ) 

That  this  is  as  much  an  expansion  in  density  as  in  small  time  can  be  appreci¬ 
ated  by  noting  that  one  obtains  nonsense  if  K^(t),K..(t)  are  expanded  in  power 
series  -  that  is  one  obtains  nonsense  if  one  attempts  to  Fourier  transform  the 
power  series  expansions. 

When  we  take  the  Fourier  transform  of  Eq .  (A6),  tne  first  term  leads  to  a 
delta  function  singularity  at  the  line  center,  u=w^ .  This  has  been  noted  in 
Baranger's  discussion  of  tne  single-pe'-tu^Der  approximation.  Moreover, 
Kf(t),K.(t)  also  contain  additional  singular  contributions  at  the  line  center. 
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The  singular-  contr i out i on s  are  connected  with  the  nonosci 1 1 atory  ter-ms  in  Eos. 
(A4)  ana  (AS),  -  .  e .  the  ter-ms  wh^c*  do  not  contain  a  sine  o'-  cosine  var-at-on 
nese  s’ngu  ar  terms  contribute  only  at  the  7'ne  center  and  are  to  oe  c-s- 
cc '•dec  -n  f~e  w-qs.  As  ciscussed  by  Bar-anger,  the  s’ ng',e-pertu'-oer  approxi¬ 
mation  is  rct  va'ic  at  the  'ire  center. 


I*  we  now  discarg  tne  nonosci 1 1 atory  terms  in  Eqs.  (A4)  and  (AS),  we  ob- 


f(t)  = - i -  e 


1 _ 


,  (A8) 

‘"Of'/S' 

where  we  note  tne  difference  in  sign  in  the  arguments  of  the  time-aeDenaent 
exponential  factors. 


With  the  above  simpl  fi cations ,  the  Fourier  transfc 


***-  -&T  I 


:an  be  performed  trivially  using 


J_  j  dt 


=  5[€aj3  +  h(u  -  ojfi)I 


(AlO) 


=  -  h(a?  -  Wfj)i  , 


(All) 
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K 


r-> 

k 


& 


k“  • 

k*.1 

.•  t." 

s.'.- 

►' 


»r*r  r»  »i  » ».7 ;  n  '.v.* '!  v  '.y,^  ■'»— ,t.*.  ** '■ "  r* 1  ■>.|nw  r»*r»T»  r* .  rr'Tr.'y 


For  $ . f(«) ,  we  ootai n 

<^jf(aj)  *  Np  E  Q{e j)  e(ek)  E  “  *  — 1 —  *  ^€ct(3  +  hoj  * 


Jk 


j  H'k' 


Ea|3 


+  E 

Jk 


2  i  <  ot '  i  V  *  j3 '  >  | 3  ^  e^fta;a'i3' 


U  e(e:)  g(ek)  U  --  ---  —  •  e  6(€a,/j, -haj  +  fia>fj)  (A12) 

^  IL  ‘  i  l»l*»  ^Ct'0'  ^ 


Ij-J'k' 


By  ma<ing  use  o'  tne  delta  functions  in  this  expression,  we  ootain 


</>if(W) 


1  r(w  -  wfj) 


7r  (efj  -  hu>)2 


(  a  :  3 ) 


where  the  naifwidtn  function  is  given  oy 


I>  -  u)fj)  =  N  E  Q(e ,)  g(ek)  E  j  <<*;V  0>  *  •  a-  5(e  +  ho;  -  ho;fi) 


Jk 

+  (w  E  e(€,)  e(€k)  E  i<a'iV!j3'>, 2  ir  6(e  ,3,  -  ho;  +  ho>fj)  (A14) 

Jk  jjj'k'  M 

This  is  identical  to  Eq.  (34)  of  tne  text  which  was  obtained  from  the  T- 
matrix  theory  of  Ref.  9.  Equation  (A13)  differs  from  Eq.  (33)  of  the  text  Dy 
a  factor  of  n  the  radiator  density. 
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appendix  b 


Miller  and  Poll  Normalization  of  Translational  States 


we  cci"3’cer  the  Quantity  in  Ea .  (37) 


=  N  E  e(€k)<k!V:k'xk'jVk>  . 


E 


<k  V  k'>  < k' !  V ! k > 


1  y  e~^ek 
n  Le 

k 


sum  = 


W  £  e‘^k 


<k,V|k'>  <k'|  Vjk> 


„  =m  /f>  the  perturber  number  density,  and 
per  per 


z  =  _ 
k 


t  E 


:he  translational  partition  sum.  With  k  /2m,  this  evaluates  to  give 


'  V  27 r/3ti2  ' 

A  particularly  convenient  set  of  basis  states  k?  for  the  present  prob¬ 


lem,  can  be  wr-'tten 


*k(F)  =  <T',k>  =  Yf  (0,  <t>)  R^(k;  r) 


where  Y  is  a  spherical  harmonic,  and  R-(k;r)  is  a  radial  wavefunction  with  the 
Miller  ard  Poll  normalization.  This  is  a  normalization  to  unity  within  a 
'aroe  but  finite  radius  Rn,  according  to 


Tr.e  asymptotic  form  of  tne  radial  wavefunction  isJ 


fy(k;  r) 


2  sin[kr  -  'frtir  +  6^(k)) 


r  —  oo 


wrere  6.,(k.)  t  s  a  orase  sn'ft. 


Next,  it  is  convenient  to  oefine  new  wavefuncti or.s  accorc’ng 

Q  ./L.  ,\  I  2  1 


R^k;  r)  =  0f(k;  0 


~om  Ec.  (B8)  the  o's  will  nave  tne  asymptotic  form  of  pure  sine  waves, 


<^(k;  r) 


r  —  oo 


sin[kr  -  V2fir  +  6^(k)] 


V>  *  ~  *>  *(VC  0  • 

■  v  Rkfm/r| ' 


(  =  ::) 


(311) 


( B12 ) 


wnere  we  have  defined 


r>  ■  &  ~~  <t>g(k;  r) 


In  Dirac  notation  we  write  this 


(313) 


|K>  =  JlT  |Rkfm>  • 


(B14) 


Now  Miller  and  Poll's  prescription  for  computing  a  "short-ranged"  quan¬ 
tity,  e.g.  a  matrix  element  of  a  multipole  interaction  whicn  vanisnes  as  r*». 


is  that  sums  over  k  can  oe  replaced  by 


E  -  E  .  E  ( 

k  ktmg  Img  0 


(E15) 
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Wnen  -q . s  (514)  and  (215)  are  inserted  into  Eq .  (B3),  we  odtain 


APPENDIX  C 


Clebsch-Gordart  Albegra 

yve  sns  ■  ■  cors'de^  t*~e  sum  (jf)  tern  as  given  oy  Ea .  (£4).  Evaluation  c* 
sun  (j.)  succeeds  in  a  similar  fashion. 


We  begin  tne  calculation  by  noting  tnat  the  three  magnetic  Quantum  num¬ 
bers  m .  .  m.  ,  .  n  do  not  cccu'-  in  the  V-matrix  elements.  Tnese  three  sums  can 
oe  oer^o^med  immediately,  ana  simdy  give 


E  E 

rrvrTij,  m 


6mj,  nij# 


•  (if  1  mfm !  if  1  jjmjXij  1  mf  ,m  j  jf  1  jjm  j() 


m 


f  * 


(Cl) 


we  use  the  aoove  Kronecher  delta  to  pe^orm  the  mfl  sum,  and  then  we  change 
some  ^ames  o*  cummy  summation  variables  to  Drovide  a  simpler  notation.  See¬ 
s''  cal  ly  ,  we  1  et 

J  ^11  J  f  ! 

'"f  i 

m^  -*  m 
m^  ,  -*•  m1 

T-'is  yields 


sum  (jf)  =  E  5(etf,  +  hw  -  e^Q  , 
jf' 

where  0  is  the  "Quantity" 

Q  =  — r  1  —  E  E  E  \v\*  . 

'z,f  +  ’  M,M'  mm' 


(C2) 


(C3) 


V  =  <jfmf;  JM;  Rkfm I V •^f'rrif '* 


C-) 


~he  oasic  strategy  for  reducing  at1  tne  multipo’e  interactions  to  a  c:r- 

,  24 

orient  form  is  to  use  the  addition  theorem  for  sphe^'ca  narmcn'c. 


Yf2m2(0-  & 


'  (2fi  +  1  )(2f 2  +  1)  1  V2 
4ir  (2t3  +  1) 


.  E 


Yf3m3^- 


•  (W2m,m2 |f ,^3m3)  .  (f,f2  00;  ftf2^30)  . 


(Oil) 


n  the  dipoie-d’Dole  case,  the  oroduct  Y.»Y,,  via  the  addition  tneorem,  pro- 

l  1 

uces  a  Y_  (which  cancels  the  f’rst  term  in  Eg.  (CIO))  and  a  Y^.  The  result  '$ 


V  =  - 


3  V4T  (1100|1120) 


.  E  E 

mrmr'  m3 


(l1mrmr,ill2m3)Y2  (0,  *>)  (#i,  •  em  )  (^2  •  em  ,)  (C12) 


New  tne  Wi gner-Eckhart  theorem  (Eq.  (20))  can  conveniently  be  put  in  an 
loerato'-  form 


m  =  Li  . 

m,  *0,±  1 


( C 13 ) 


^here,  'f>om  Ea .  (20),  the  matrix  element  of  tne 


u(m, )  are  gi ven  by 


(C14) 


- rom  Eq.  (C13)  we  then  note 


f1  *  €mr  U  ^  emi  •  em  . 

mi  =0.  ±  1  r 


=  E  fi{ m,)  (-1)  r  5 


mi,  -mf 


M  •  e_  »  (-1)  mHV  • 


(Cl  5) 


The  di po  i  e-di pole  interaction  then  assumes  its  final  form 


V  =  - 


3  v47T  (1100  : 1120) 


E  E  (11mrmr,!112m3)Y2m3(0,  <p) 

mrmr*  m3 


mr  +  rtir/ 

v  1  r 


(C16) 


In  this  form  tne  matrix  element  is  easily  evaluated;  operates  De¬ 

tween  R  -  -  states,  y,(-m  )  ODerates  Detween  j^m^  raciator  states,  anc 
k  1  m  1  r  1 

p0(-m  ,)  operates  between  JM-  perturber  states.  The  complete  matrix  e  ement 
u  r 

V  in  Eq.  (C4)  can  then  be  written 


V  =  -  3  v'4t 


(110011120) 


n  tn  m. +  mr, 

L,  u  (-1)  r  (11  mrmr,|112  m3) 

mrmr,  m3 

’  <YfmiY2m3^*  ^|Yf'm'> 

•  <  jfmf  |  /i,(-mr)  I  jf '^f '  >  <  J  M !  nd-mr,)  j  J '  m  '  > 


(C17) 


Matrix  elements  of  the  y  operator  are  obtained,  in  terms  of  reduced  matrix  el¬ 
ements,  using  Eq .  (C14).  Matrix  elements  of  the  Y  operator  are  ootainec  from 

.  ,  16 
the  formula 

<YfmlYf,, n,3<«.  *»l 


=  [~  ~4T(2f 1 (f 3* 'mjm ' ■  *  {M’OQ.lit'tO) 


(C18) 
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Tne  aoove  formula  obtains  by  combining  the  addition  theorem.  Eg.  (Cii).  wito 
tne  crthcgcna 1  i  ty  '•elat-on  for  soherica1  harmonics 


coon  souar-ng  tne  matrix  element  in  £q  .  ( C 1 7 )  and  inserting  ;nto  Eq . 
(C3),  tne  sums  ove"  the  m-indices  can  be  carried  out  using  one  o*  tne  *ur-ca- 
mentai  orthogona 1 i ty  re'at’ons  for  Clebsch-Goraan  coefficients,  namely 


E 

mtm2 


(jij2mim2ij1j2j'm')(j1j2m,m2:j,j2jm)  =  5j,j<$m,m' 


o<"  for  j '=;  and  m'=m 


E 

m  i  m2 


(iiiamim2  =  1 


( C2C) 


(-21) 


In  carrying  out  this  exercise,  one  has  to  make  use  'of  symmetry  proDerties  of 
CG  coefficients. *  * • 1 3 

Tne  final  results  for  the  dipole-dipole  case  is 

2 

Q  =  —r~  (2f  +  1)(2J  +  1)(f200  |f2f'0)2 

#l<Rkfl  lRk'f'=>  2 

,  ( C22) 

•  I  <JfMMi!!jr>l2  !  <J!  'nil  !J'>  i* 

liDon  inserting  Eq.  (C22)  into  Eq .  (C2),  we  arrive  at  the  result  for 
sum  (jf)  quoted  in  Eq .  (44)  of  the  text. 


Cl po 1 e-Qua drupole  Case 


The  basic  strategy  of  using  the  addition  theorem  Hs  the  same  as  before 
but  a  good  deal  more  comqlicated  to  carry  out.  We  present  only  a  few  of  the 


7?  *\>  r.*  T'-*  -r  <-■  *■ 


rne  ci pel e-quaarupol e  interaction  can  be  written 


{2  f  •  Q  •  jt  -  5{r  •  /Z)(r  •  Q  .  f)] 


c::  i 


In  this  expression  we  snail  taxe  y  to  be  tne  aipcie  moment  o'  tne  '•aerator  arc 
Q  to  be  tne  quadruooie  moment  tenser  of  the  perterber.  The  '■everse  success 
also  occurs  ana  nas  peen  inciudec  in  our  calculations. 


‘ne  traceiess  Quacrupole  moment  tensor  is  Given  by 


Q  =  E  e)(f|r]  -  -L  i). 


( C24 ) 


wnere  I  is  tne  unit  tensor. 


in  reducing  Eg.  (C23),  we  use  Eq .  (C6)  for  r,  tne  operator  form  of  tne 

Wigne—Ecknart  theorem  embodied  in  Eqs.  ( C 1 3 )  througn  (CIS)  for  J,  anc  fi- 

« 

naily,  a  particularly  useful  tensor  Wigner-Eckhart  theorem  for  C.  namely 


m3  =  0, 1 1,±2 


( C25) 


m,  *0,  ±  1  m2  *  0,  ±  1 


(1 1m,m2i  1 12m3)  (C2  6) 


wnere  the  matrix  elements  of  the  operator  G0D("Mp  are  given  by 


<  JiMi ;  QOQ(-ni3)  I  >  ■  <  J|  i !  Q |  j  J2 >  (J2  2  -  M2  -  rrij  |  J2  2  Ji  -  Mi) 


( C27) 


and  where  'J,  Q  is  a  reduced  matrix  element 


By  -"eDeated  use  of  the  addition  theorem  to  Eq.  (C23).  the  di pole-quadru- 

25 


sole  interaction  can  oe  reduced  to  tne  form 

V  =  - 


4t  15 


1 

r4 


3  7 

SEE  (-i)mrM(-m) 

mr  =  0, ±  1  m3  =  0.±l,±2  m"  =  Ot±l 


•  (-1)m3Qop(m3)(12mrm3:i2  3m")  •  Y3m„(0,  <p)  (ZZS) 

”ne  matrix  elements  of  V  above  can  now  be  obtained  using  Eos.  ( C  i  ) , 
(CIS),  and  (C27).  Upon  soua'-ing  the  matrix  element,  inserting  into  Eq .  (23), 
and  carrying  out  the  m-summations  using  Eas.  (C20)  and  (C21),  we  find  as  the 
(i'na'  '•esu’t  for  tne  ci  do  1  e-auaoruoc '  e  case 


Q  «  (~~f  0  +  1)(2J  +  1) 

.  (tf300!f3f'0)2 

•  <Rkr  ~r  '2 

•  !  <J| ! Q !  IJ'> I2 


~n i s  can  be  readiiy  be  comoared  to  Eq.  (C22)  for  the  di poi e-di do! e  case.  The 
differences  are 


( C30a) 


1 


r3 


( C3Cb) 


(f200  f2f'0)J  — ►  (f300  ftf'O)2 


( C30c ) 


<J|  \fi\ ! J' >  i 2 


I  <J| ! Of  |J'>  |2 


( C30d) 


I'  ^ . 


i 

■? 

§ 
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